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The generic properties of the emission of coherent radiation from a moving charge distribution are discussed. The general structure 
of the charge and current distributions in an extensive air shower are derived. These are subsequently used to develop a very intuitive 
picture for the properties of the emitted radio pulse. Using this picture can be seen that the structure of the pulse is a direct reflection 
of the shower profile. At higher frequencies the emission is suppressed because the wavelength is shorter than the important length 
scale in the shower. It is shown that radio emission can be used to distinguish proton and iron induced air showers. 
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1. Introduction 

In recent years the field of radio detection of cosmic ray air 
showers has reached a mature stage which is shown by the many 
contributions to this meeting on the topic of radio emission. 
Consensus is approaching on the theoretical description |[T) of 
the emission process and a detailed quantitative understanding 
of the experimental results (H O is close although some chal- 
lenges remain [4|. Installation is in progress of extensive arrays 
of radio detectors at the Pierre Auger Observatory [5 |, at LO- 
FAR [6], and at the South pole Q. In this work the importance 
of coherent radio emission from air showers will be stressed 
where the most important emission mechanisms were already 
investigated in the earliest works on this subject ll8ll9l [T0l[TTI . 
namely Cherenkov and geo-magnetic radiation. A complete 
historical review is given in Ref. ifTZl . 

Coherent emission occurs when the emitting charges are con- 
fined to distances A which are smaller than the wavelength A. 
Since the source is 'viewed' with a resolution of the wave- 
length, a fine sub-structure in the source will not affect the 
emission process. At a much shorter wavelength, A «: A, 
the different parts of the fine structure in the charge distribu- 
tion will contribute to the emission process with varying phases 
which as often result in constructive as in destructive interfer- 
ence with as result that the net emission probability is strongly 
suppressed. In an EAS initiated by a cosmic ray of 10 18 eV the 
number of charged particles at the shower maximum is of the 
order of = 10 6 and coherent radiation, where the intensity 
is proportional to A^ 2 , is far more intense that incoherent radia- 
tion where the intensity is proportional to N. Only at high fre- 
quencies, where the coherent process is suppressed because the 
wave length is much smaller that the relevant size of the emit- 
ting charge distribution, the incoherent process can contribute. 
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For coherent emission thus only the macroscopic structures in 
the EAS contribute. 

At high frequency the coherence condition A < A will no 
longer be satisfied implying a cut-off of the coherent response. 
Some of the length scales that are important for EAS emis- 
sion [ 13 1 are i) the pancake thickness; ii) the length of the EAS 
projected along the line of sight; iii) the lateral distribution of 
the charges in the EAS. The high-frequency cut-off is reflected 
in the time between the start and the zero crossing of the pulse. 
The important aspect of the macroscopic model is that it clearly 
shows which are the important lengths that determine the emis- 
sion process. 

Model independent conclusions can also be drawn at large 
wave-length since the intensity of the emitted radiation of a sys- 
tem of charges decreases linearly with increasing wave-length 
when considerably larger than the the size of the emitting body. 
An EAS, independent how enormous the event may be, always 
exists for a limited time and occurs in a limited part of the at- 
mosphere and currents are confined to a limited region of space- 
time. As a consequence the intensity of the emitted radio waves 
should vanish linearly in the limit of infinite wave length or zero 
frequency. The time integral of the emitted pulse should vanish 
implying equally large positive and negative amplitudes. The 
simplest structure of a pulse is thus bi-polar and unipolar pulses 
are un physical lfT3l . 

For a correct description it is important to have a consistent 
description of the charges and their motion. Consistent in this 
respect means that all charges are accounted for and thus charge 
conservation holds. Charges may move and be suddenly accel- 
erated (through a collision) but no net charge can be created. 
For example, if an electron is accelerated in a Compton scat- 
tering process, a positively charged ion should remain behind. 
Once the consistent description is obtained of the distribution of 
the charges and their velocities, resulting in a four-current den- 
sity /*(/, f), the emitted radiation is straightforward to calculate 
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by applying the Maxwell equations with the current density as 
input. As mentioned, for coherent emission only densities av- 
eraged over an appropriately chosen length scale matter where 
one should be careful to account for all charges and currents 
in the system to obey charge conservation. Due to the neces- 
sary averaging over the path of many individual electrons a very 
simple picture emerges for the emission process which lies at 
the basis of the Macroscopic GeoMagnetic Radiation (MGMR) 
model (HE). 

In Section [2] the essential current distributions for a generic 
EAS are derived which form the basis for the (semi) analytic 
predictions for the structure of the radio pulse in the MGMR 
model lfl4l [T31 . As an application we show, using a hybrid ap- 
proach where the parameters of the MGMR model are extracted 
from a Monte-Carlo simulation of the air-shower, that the radio 
emission can be used to distinguish proton and iron induced 
showers J4). 



2. The Macroscopic Model for Geo-Magnetic Radiation 

The pancake at the front of an EAS consists out of a plasma 
with large amounts of electrons, positrons and other particles 
moving towards the surface of the Earth with a velocity almost 
equal to the light velocity. In the center of the shower con- 
stantly electron-positron pairs are created, mainly by the ener- 
getic photons in the shower core, to form the pancake, depicted 
by the broad band in Fig.[TJ The shower front, as it is driven by 
photons, moves towards the Earth with the light velocity while 
the produced particles trail some distance behind the front and 
are responsible for the finite thickness of the pancake. A net 
electric current is induced in the pancake because of the Earth's 
magnetic field induces a Lorentz force pulling the electrons 
and positrons in opposite direction. Because of the constant 
interactions with the air molecules the leptons (electrons and 
positrons) reach a constant drift velocity Vd as corroborated by 
Monte-Carlo shower calculations |4|. In Fig. [T]the drifting elec- 
trons and positrons are denoted by the chequered arrows. The 
motion of the electrons and the positrons thus contribute coher- 
ently to a net electric current density in the direction indicated 
by the red arrow. Due to the constant interaction with the ambi- 
ent air molecules the leptons will loose energy and trail further 
behind the shower front and become non-relativistic. In the fig- 
ure these leptons are denoted as 'stopped'. On average the elec- 
trons and positrons are separated by a distance D$- The fact 
that this happens can also be seen as a consequence of charge 
conservation, the electric current must induce a displacement of 
net charge from one side of the shower core to the other. One 
consequence of the induced current density may not be directly 
obvious. As the electrons and positrons move in opposite direc- 
tions there is already a net displacement of charges that move 
with relativistic velocities. Their displacement is about half of 
the stopped charges, D M ~ As /2. This completes the qualita- 
tive picture of the charges and currents that are induced by the 
geomagnetic field. For the following discussion it is instructive 
to distinguish the different components which must be present 
in any realistic shower simulation. 
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Figure 1: [color online] Schematic description of the current densities in an 
Extensive Air Shower 1141 . 



In addition to the current and charge densities of geomagnetic 
origin there is also an induced charge that is independent of the 
magnetic field. Through the process of Compton scattering on 
electrons in air molecules and electron knock-out reactions by 
the relativistic electrons a net excess of relativistic electrons is 
created in the pancake. Simulations indicate that this excess 
can be large, about 30% of the total lepton density Q and this 
is indicated in Fig. [T]by the star in the shower core. The net 
negative charge excess implies that there must also be a positive 
charge in the system, which is formed by positively charged air 
molecules. Even though these positive charges are at rest, they 
contribute to the radiation field since their number increases as 
the shower develops, but are not marked in Fig.[T]for simplicity. 

To obtain better insight in the emission mechanism we will 
use the picture of the currents described in Fig. [T[ realizing 
that all charges are concentrated very near the shower axis. As 
is discussed in detail in Ref. IIT41 [151 this results in a picture 
where the geomagnetic electric current becomes point-like at 
the shower axis and the charge separation of the relativistic lep- 
tons is taken into account as an electric dipole moving at a rela- 
tivistic speed. The stopped particles are similarly accounted for 
by an electric dipole that is at rest in the Earth frame but con- 
stantly increasing in magnitude as the shower front proceeds to 
the surface of the Earth. As shown lfl4l[T5ll . the induced electric 
current and the net charge excess give the dominant contribu- 
tion to the emitted radiation. 

To obtain an analytic expression for the pulse due to the ge- 
omagnetic current the thickness of the pancake is assumed to 
be small 1 14 1 and we adopt a simple geometry with a verti- 
cal shower and a horizontal magnetic field. The number of 
electrons and positrons in the shower at an height z = -ct r is 
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parameterized as N(z) = N e f t {t r ) in terms of the normalized 
shower profile, f r (t r ), where N e is the number of electrons in 
the shower at the maximum. The induced geomagnetic current 
in the i-direction is 



j(x,y,z,t) = (v d q}eN e f,(t r ) 



(1) 



In the simple picture the drift velocity is assumed to be inde- 
pendent on the height in the atmosphere. In the limit where the 
shower moves with the light velocity and the index of refrac- 
tion of air equals unity, the retarded time t r can be expressed in 
terms of the observer time t as 



ct r « — - — 



let 



(2) 



The observer is at a distance d from the core and t is time after 
the shower hits the surface of the Earth. Eq. Q shows that the 
early part of the received pulse is emitted at large (and negative) 
retarded times and thus large heights while the late part of the 
pulse is emitted when the EAS was already close to the round. 

The only non-vanishing component of the vector potential is 
in the direction of the electric current, 



A x (t , d) = J 



ftttr) 
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where J = (vdq)N e e/4neoc is a constant depending on the en- 
ergy of the cosmic ray and T> is the retarded distance. The 
electric field is proportional to the time derivative of A, giving 



c 2 t 2 4 d 
E x (t,d)*J—^—[t r f t (t r )], 



(4) 



d 4 dt r 

using Eq. |2]). From this simple expression some important ob- 
servations can be made which also hold for realistic calculations 
to a surprisingly accuracy. The first is that radiation is emit- 
ted due to the variation of the shower profile with height. The 
zero crossing of the electric field at the observer position corre- 
sponds to the maximum in t r f,(t r ) and thus the emission from 
the shower at an height exceeding that of the shower maximum. 
The dominant part of the pulse is thus emitted at heights well 
above the shower maximum which implies that the radio signal 
bears information on the early stages of the EAS development. 
A second observation is that, because the time dependence of 
the pulse is expressed in terms of the shower time t r Eq. it 
can be seen that at twice the distance from the point of impact 
the pulse is four times as wide while the amplitude is decreased 
by a factor 2 4 . It can easily be seen that this simple picture is not 
valid at small distances to the shower core where the pancake 
thickness becomes the determining length scale. 

A prime example showing the strength of the macroscopic 
approach is the introduction of the drift velocity. The motion 
of any given electron (or positron) in the shower is very com- 
plicated. It is suddenly accelerated (or created) and moves at 
relativistic velocities in the Earth magnetic field where it expe- 
riences a constant acceleration due to the action of the Lorentz 
force. Through soft processes, mainly Thompson scattering, 
it is constantly loosing energy and decelerating. It may be 
strongly deflected in hard scattering processes. With all acceler- 
ation processes radiation will be emitted. For coherent emission 
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Figure 2: [color online] Drift velocities as function of shower time t r for a few 
proton (blue curves) and iron (red curves) initiated showers. 



the motion of all electrons in a certain phase-space element has 
to be averaged. At any time some of the electrons in a volume 
will be accelerated by the Lorentz force while others are decel- 
erated by the action of the air. The net result is that on average 
the electrons move with a constant drift velocity in one direc- 
tion and the positrons in the other. This picture is supported 
by Monte-Carlo calculations showing that the drift velocity de- 
pends on air density as shown in Fig. [2] As a result of reaching 
an equilibrium with an (almost) constant drift velocity neither 
synchrotron emission (acceleration) nor bremsstrahlung (decel- 
eration) is of particular importance and only the time variation 
of the induced macroscopic electric current causes the EAS to 
emit radiation. 

In the complete MGMR calculations a 4-current, j^if, f), is 
constructed where the zeroth component corresponds to the 
charge and the other components to the electric current. The 
electric field at the observer position d and time t can be ex- 
pressed as a derivative of the Lienard-Wiechert potential which 
in turn can be calculated directly from the 4-current density fol- 
lowing any standard textbook on Electromagnetism |[T6l 
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ignoring the lateral extension of the shower profile and where 
the integration runs over the pancake thickness h. We use the 
common notation where h is a unit vector pointing from the 
source to the observer and R is the distance, both evaluated at 
the retarded time. 

3. Distance dependence of the radio pulse & composition 

Up to this point we have focussed on obtaining a qualitative 
description of the current densities in an EAS. To make quanti- 
tative comparisons with data a hybrid approach has been devel- 
oped where the important model parameters are obtained from 
Monte-Carlo simulations |4|. 
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These simulations are done using the cascade mode of the 
CONEX shower simulation program (where further details can 
be found in Ref. JU) modified to include the deviation of 
charged particles in the Earth's magnetic field fl5l . An analysis 
tool has been written to give the full three-dimensional and tim- 
ing information of the currents and particle distributions in his- 
tograms. More details of the complete shower-simulation pack- 
age including analysis tools, called CONEX-MC-GEO, will be 
discussed in a future publication. 
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Figure 3: [color online] Total number of electrons and positrons as function 
of shower depth for 40 proton (blue curves) and 40 iron (red curves) initiated 
showers for E = 10 18 eV. 

As an example CONEX-MC-GEO, is used to study the dif- 
ferences of the Lateral Distribution Function (LDF) for two 
types of showers, proton- and iron-induced. Simulations are 
performed for 40 proton- and 40 iron-induced showers. The 
shower profiles show much larger shower-to-shower fluctua- 
tions for proton- than for iron-induced showers, see Fig. [3] As 
also is expected, the penetration depth, X mnx , is larger for proton 
than for iron. The simulations are done for the site of the Pierre 
Auger Observatory where ground level is at 870 g cirT 2 . The 
fraction of excess electrons not very dependent on initial energy 
or composition of the cosmic ray but is slowly increasing with 
shower depth. Near the shower maximum the charge excess 
is close to 23% of the total number of electrons and positrons 
while at ground level it has increased to 25%. As a measure of 
the pancake thickness we have calculated the mean distance of 
the electrons behind the shower front, L = (h), the mean pan- 
cake thickness parameter. This shows that the mean pancake 
thickness parameter is almost independent of shower height or 
energy. It depends on the shower type, varying from L — 3.9 m 
for proton-induced air showers to L = 4.3 m for iron. 

The extracted parameters are used in MGMR calculations. 
The maximum of the field strength as function of distance is 
shown in Fig.|4]for an observer placed on the x-axis where the 
charge-excess and geo-magnetic fields interfere constructively. 
The slope of the LDF for iron-induced showers is less steep 
than that for protons. Close to the shower core the pulse struc- 
ture is determined by the pancake thickness which is the reason 
for the minor differences between the two shower types at small 
distances. Only at large distances the pulse structure is deter- 
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Figure 4: [color online] The LDFs are shown for proton (blue diamonds) and 
iron (red crosses) induced showers at an energy of 10 eV, for an observer po- 
sition at a position where the charge excess and the geo-magnetic contribution 
interfere constructively. 

In a real experiment the pulse will be filtered in frequency. 
For this reason we calculated a ratio that is closer to the data, 
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P(50,/>25)/P(300,/>25) 



(6) 



where P(d, f > 25) is the power in the pulse for frequencies 
in excess of 25 MHz at a distance dm from the shower core. 
The value of ^5o/ 300 is calculated for each of the simulated 
proton and iron induced showers separately. Shower-to-shower 
fluctuations result from the fact that for each of the simulated 
showers a different profile is extracted. The results are dis- 
played is Fig. [5] as a histogram. The averages for iron- and 
proton-induced showers show a large difference, however due 
to shower-to-shower fluctuations some of the proton showers 
result in similar values as iron, directly reflecting their shower 
profiles. At an energy of 10 17 eV the average value of ^5o/ 300 
differs more. This is a direct consequence of the fact that the 
important parameter responsible for the observed effect, X max , 
differs more. 

Due to interference of the charge excess and the geomagnetic 
contributions the LDF, and thus also ^5o/ 300 , depends on the ori- 
entation of the observer with respect to the shower axis. It has 
been checked that the qualitative dependence on shower type 
remains as shown in Fig. [5] this because the the charge-excess 
fraction is similar for the two types of showers. 



4. Summary and Conclusions 

We have shown that the Macroscopic GeoMagnetic Radia- 
tion model for coherent radio emission from extensive air show- 
ers offers a very realistic picture for the emitted pulse. It is sim- 
ple enough to give hands-on insight in the emission mechanism 
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Figure 5: S^/ 3 qqi see text, for the 40 simulated proton and iron showers at an 
energy of £ = 10 18 eV. 

which is of help for the interpretation of the experimental data. 
Depending on the distance from the shower core there are dif- 
ferent length scales that determine the structure of the coherent 
radio pulse. Close to the shower core the pulse duration is de- 
termined by the pancake thickness while at larger distances the 
shower profile is essential. This offers the possibility to use the 
LDF of the radio pulse to distinguish iron and proton induced 
showers even though the effects of shower-to-shower fluctua- 
tions for protons can be large. 

In this presentation little attention has been given to the pos- 
sibilities to use the polarization of the emitted radiation to dis- 
entangle the charge excess and geomagnetic components (4). 
Preliminary experimental results look very promising [3 1. Also 
the effect of the lateral distribution of the particles in the shower 
and that of a finite index of refraction of air [ 15 1 are being in- 
vestigated. 

5. Acknowledgment 

This work is part of the research program of the 'Sticht- 
ing voor Fundamenteel Onderzoek der Materie (FOM)', which 
is financially supported by the 'Nederlandse Organisatie voor 
Wetenschappelijk Onderzoek (NWO)'. 

References 

[1] T. Huege, M. Ludwig, O. Scholten, K.D. de Vries, these proceedings 
(2010), |aFXrv7l 009.0346 i2 [astro-ph.HE], 

[2] M. Ludwig, T. Huege, these proceedings (2010); Astropart. Phys. submit- 
ted (2010). 

[3] H. Schoorlemmer [Pierre Auger Collaboration], these proceedings 

(2010) ; D. Fraenkel [Pierre Auger Collaboration], Nucl. Instr. Meth. A 

(20 1 1 ) , Proceedings of ARENA 20 1 0. 

[4] K. D. de Vries, O. Scholten, K. Werner, these proceedings (2010); K. D. 

de Vries, A. M. van den Berg, O. Scholten, K. Werner, Astropart. Phys. 

in press (2010), |arXrv:1008.3308| 
[5] S. Fliesscher [Pierre Auger Collaboration], these proceedings (2010). 

5 




